Thermal conductivity of a granular metal 
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Using the Kubo formula approach, we study the effect of electron interaction on thermal transport 
in the vicinity of a metal-insulator transition, with a granular metal as our model. For small values 
of dimensionless intergrain tunneling conductance, g <C 1, we find that the thermal conductivity 
surprisingly shows a phonon-like algebraic decrease, k(T) ~ g'^T'^ / even though the electrical 
conductivity obeys an Arrhenius law, a{T) ~ ge~^''^'^; therefore the Wiedemann-Franz (WF) law 
is seriously violated. We explicitly show that this violation arises from non-magnetic bosonic exci- 
tations of low energy that transport heat but not charge. At large values of intergrain tunneling, 
we find it plausible that the WF law weakly deviates from the free-electron theory due to potential 
fluctuations. Implications for experiment are discussed. 



At low temperatures, transport in many metals is dom- 
inated by elastic scattering of electrons. In the absence of 
phase transitions (such as superconductivity) and collec- 
tive modes, the Wiedemann-Franz (WF) law [J relates 
the electrical (tr) and thermal (k) conductivities through 

a universal Lorenz number, Lq = Ta(T) ^^^^^ = 3^2° 
is Sommerfeld's value for the Lorenz number. Devia- 
tions from the Wiedemann-Franz (WF) law beyond those 
due to phonons and inelastic processes such as electron- 
phonon scattering indicate a presence of collective modes 
that discriminate between charge and energy flow. Such 
deviations are expected to get more pronounced as one 
approaches a metal-insulator transition. 

The issue of excess electronic thermal conductivity in 
disordered metals with interacting electrons has been a 
topic of much theoretical debate in recent years. The 
WF law has been afSrmed for any strength of scattering 
0, disorder Q, or magnetic field 0. In Refs. 0,01 1 the 
WF law was shown to hold on account of a Ward iden- 
tity. Subsequently there have been many studies includ- 
ing the effect of electron interaction, some of which have 
suggested deviations from the WF law [TllllglllollTllll^ . 

In this paper we study the effect of Coulomb blockade 
on the thermal conductivity of a nonmagnetic granular 
metal. Granular metals are a novel playground, experi- 
mentally and theoretically, for studying electron interac- 
tion effects in the vicinity of a metal-insulator transition. 
At not too low temperatures transport in a granu- 
lar metal is due to a competition of incoherent intergrain 
tunneling and Coulomb blockade [13, Q. The result- 
ing model captures the physics of coupling of the contin- 
uum of electron states in the metallic grains with electro- 
magnetic fluctuations. Granular metals differ from Hub- 
bard models in having a near-continuum of excitations 
in each metallic grain. Their study is relevant beyond 
artificially-prepared granular metals 0, 0| . For exam- 
ple, in more insulating regimes, these materials show a 
similar electrical conductivity to disordered semiconduc- 
tors in the VRH regime. In more conducting regimes. 



granular metal physics appears to describe transport in 
underdoped cuprate superconductors llTtllSj where there 
is electronic phase separation into metallic and insulating 
regions (l9| . 

We summarize our main results. We obtain for the first 
time an explicit Kubo formula for the thermal conductiv- 
ity of the AES model through a double functional differ- 
entiation of the partition function with respect to an ap- 
propriate source field jl^ • For small values of the dimen- 
sionless intergrain conductance, g <^ \, we identify low- 
energy non-magnetic bosonic excitations that contribute 
to heat transport but not charge transport. Conse- 
quently, as our calculations show, the thermal conductiv- 
ity decreases only algebraically, k(T) ~ g^T^/E^, quite 
unlike the electrical conductivity of a regular granular 
array, which obeys an Arrhenius law, cr(T) ~ ge~^''^'^ . 
These bosonic modes are the nearly- continuous particle- 
hole excitations in the metallic grains and have no coun- 
terpart in a half-filled Hubbard model with U ^ Ec where 
particle-hole excitations are gapped. Since the electrical 
conductivity is insensitive to these neutral excitations, 
the algebraic thermal conductivity we calculate distin- 
guishes between a granular metal and a Hubbard model. 
This finding also reminds us of the excess heat trans- 
port through excitons proposed recently H^l in the 
context of disordered semiconductors. For strong tunnel- 
ing, g 3> 1, we find it plausible that the Lorenz number 
weakly deviates from the Sommerfeld value Lq due to 
Coulomb interaction. 

Consider the following microscopic Hamiltonian for a 
granular metal array ignoring spin interactions. 



xA 



^ = E V4exA^xA + ^c5I(Qx-Oo 



(1 + h^) + 



\ E (tii'^lAV'x'A'+h.C.)(l + /!.), (1) 



{xx'),AA' 



where x labels the grains and A labels the states within 
each grain, ^xA = fxA-/^- Qx = Y^x '/'IaV-'xA is the charge 



2 



on grain x, Ec is the grain charging energy, and t^^, is 
the matrix element for nearest-neighbor hopping. Small 
differences in temperature of different grains are modeled 
through the introduction of "gravitational" potentials /ix 
[201] that couple to the grain electron energy, li 13 — \/T 
is the average inverse temperature, then hy^ — (/3x — P) 
is the small excess of (3 on the x*'' grain. 

The strong Coulomb interaction in Eq. makes it dif- 
ficult to calculate correlation functions directly using the 
electron operators. We therefore decouple the interaction 
to obtain an action for a Hubbard- Stratonovich field, V , 
representing the electrostatic potential on each grain, 



respect to /irx- Dropping any constant terms, and using 
the above Green functions, we have. 



^ \t\A[ G{T„r)x 

I J 



, ■_, IJTI 

x'— x±a 

X G(t, Ti) C0s(.^^xx' - (ySrixx') + ^ CxaGxa(ti, t) 



X Gxx{t, T2)G^'y{T2,Ti) exp[i(93r2xx' - (/^Tixx' 



-) + 



1 f 

xx';AA' 

+ V / dT^ija,-fy,x + exA(i + /irx))^xA, (2) 

xA 

where r is imaginary time, and the h-x_ sources are now 
T— dependent. We then make the gauge transformation 
"^xA — * e~"^^="!/'xA and 14; = idr'fTx- This ehminates 
the tp^V^p term, at the cost of replacing t by a gauge- 
dependent hopping amphtude, i^^^, = t^^,e^'^^^^' , where 
V'rxx' = frx — frx'- Thc Matsubara fields lyS^x satisfy 
bosonic boundary conditions, (pf}^ — 27rfcx + </'0x, where 
/cx G Z is the winding number at site x. 

Next we integrate out the fermions and expand the 
electron determinant to 0{ht^), 



where we assumed that the bare tunneling element is a 
constant. To obtain an expression for the energy cur- 
rent J7-XX' from grain to a neighboring grain x', we 
employ the continuity equation, drE^.^ — iJ2x''^rxx'- 
We will need the equation of motion for (^^tx from Eq. Q 
(with sources set to zero) and the Green functions de- 
fined in Eq.Q. The symmetrized energy current, j'l^^, — 

2('^TXx' Jtx'x)j is 



Jrxx' = ^^gT 

X sin(y)^xx' 
1 



'-[dr'Prx + dr^Tx')- 



1 



dr 



sin[7rr(T — t')] 



sin[7rr(T — t')] 



— r sin((^.rxx' - ^Pt'xx'), 



(6) 



(drifrxf 



+ tr[GxA(T, T)/lrxCxA] 



where we have introduced the dimensionless intergrain 
tunneling conductance g — 2Triy'^ {ep)\t\'^ and used 
G(t,t') — . ^"^f'''^ ,K . To obtain the Kubo formula for 
thermal conductivity, we need to calculate the non-local 
specific heat through a double differentiation of the par- 
tition function with respect to hrx, and then extract the 
thermal conductivity using the continuity equation. 

The same result may be obtained more easily by in- 
troducing a source field frx that can generate the energy 

+tr[GxA(ri,r,)/i.,xexAGxA(r2,r3)f^3^xx'Gx'A'(r3,r0f^:iJ,(3)™.t' In Eq.® we expand to linear 

order m hrx and use the held jrx instead 01 hrx- Next 
we integrate out the conduction electrons and expand to 
second order in tunneling. The resulting action is 



X 

+ ^tr[GxA(Ti, r2)i:^^'',,^,Gx'A'(T2, ri)t 



+A'A 1 
rix'xJ 



+ xtr [GxA (t 1 , T2)t^2xx' Gx' A' (t-2 , Ti )f^jx'x (^Irsx 



^T2X' ) 

wA'A 



where "tr" means a sum over all indices, and Gx.x{t,t') 
is the electron Green function 



— (r — r') 



GxA(r,r') = T5:— - 



xA 



(4) 



Here — 2TrT{n + 1/2) and Qr is the unit step func- 
tion. We will also need the momentum-summed Green 
function, G{t,t') — X^a GxaI^, t'), which is indepen- 
dent of the grain label. If the temperature is much 
larger than the Thouless energy for intergrain diffusion, 
T » |tpi^(eF), and if the grains are much larger than the 
Fermi wavelength of the metal, it suffices to expand the 
electron determinant to second order in tunneling 23] . 

The (time-dependent) energy of the x*'* grain, Erx, 
is obtained by differentiating the action in Eq.® with 



+ 2 tr [GxA (ri , )t>;^^^, G^,y {r^ , n Kx'x] , (7) 

t^rxx' = i^rxx' [l + */rxx' ((CxA + Cx'X')/2 + ^i™))] , (8) 



where /,xx' = frx-frx' and l/J™? = {Vrx + Vrx')/2. The 
energy current in Eq.© can be obtained from Eq.Q us- 

ing j).^^,[ip,f] = SS[ip, f]/5frxx'- This is similar to the 
derivation of the electric current, with electric charge re- 
placed by electronic energy. The thermal conductivity k 
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may be obtained by further differentiation, 



1 

UJ Jn 



Dip e 



E 



■{E) 



'^irfi,x+a[yJ] 

^ frix^x-l-a. 



'r2x',x'+a J y 



Tix' ,x'+a 



,(9) 



,(10) 



where U!~^ ~ lu + iS. The first term, which has come to 
be known as the "diamagnetic" contribution is a 

local term in the sense that it is an average on a single 
bond. The second term in Ea.l|l()|l is a product of two 
energy currents, and is also known as the "paramagnetic" 
contribution k^^^. The diamagnetic contribution may be 
obtained by differentiating Eq.Q twice with respect to 
the source. 



^ Jo 

sin itTt sin 7rT(T+i-rc) "^rx.x+a "^Ox.x+a sin^ ttTt 



X (cos (<^TX,x+a - </20x,x+a)) 



(11) 



,rc — ^0 



Let us call the independent part of Ea. (|ll|l the ki- 
netic contribution, k'-'^'^\ and the ^-dependent part the 
potential contribution k^'^'^K Ea. (|ll|l is similar to the 
known expression for the diamagnetic contribution to the 
electrical conductivity obtained in Ref . [2^ . 

A. Consider first weakly coupled grains, g <C 1. Do 
a perturbation expansion in g, k = k'^^^ + k*^^^ + • • • , 
where The diamagnetic parts of the electrical 

[i^ l and thermal conductivities [Eq. (fTT] involve the bond 
correlator, 11^ = (exp i((/j^x,x+a - <y5ox,x+a)) iK") + 
n^^) + • • • . Because of the Coulomb blockade, H?^ = 
lEgi 9, e-(«i+'?2)'3£;<^-2Scr(i-5i-?2) jg exponential in r. 
So 

In contrast, n^^'' has a power-law r-dependence l25f . 
which is a result of the l/r^ interaction in the AES ac- 
tion, Eq.Q, arising from particZe-Zio/e inelastic cotunnel- 
ing: 



(1) ^ 2ngT^ ^ 2t^ 



sin (ttTt) 



It gives rise to 0{g^) terms in both the diamagnetic and 
paramagnetic contributions to a and k. We have shown 



[2 51 in a previous work that 

27rg 



(e' 



l)a^ 



r2,, — *— 40; + 



47r 



the power-law terms in the electrical conductivity can- 
cel [2^, leaving an Arrhenius law a = cr'^'^'^-' e~^=/"^. 
Now consider the 0{g^) terms in the thermal conductiv- 
ity. The dominant contributions to both the diamagnetic 
and paramagnetic terms arise from charge-neutral config- 
urations {qi =92 = 0), and, unlike the 0(5) terms, are 
not Arrhenius-suppressed: 



\ / ^ r 



(13) 



where 7^ = T'^{s\i\{-KTT))-'^dl{sa\{-KTTy^ . These two 
terms have the same sign and do not cancel^ 



12^3^2_, (g^T^ 
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(14) 



and so the thermal conductivity is only algebraically 
small, 



^(K.2) ^ (l67rVl5)(.g^TV£^c^). 



(15) 



This is the central result of our paper. In the argument 
above, the algebraic temperature dependence is an even- 
tual result of the particle-hole inelastic cotunneling pro- 
cesses described by the AES action. 

The potential energy contributions k^^^^^ 
and K^^'^^ which involve averages such as 



(K-x,x+a^0x,x+aCOs((^r 



- <;50x,x+a)), turn out 



to be proportional to (91 — 92) • This immedi- 
ately rules out neutral processes, such as inelastic 
cotunneling, for which 91 = (72 = 0. The 0{g'^) con- 
tribution to the above average can be shown to be 

EcEg,,gM - g2)^e-(9?+«')/3S.-2£;.r(l-5i-,,)^ ^^^^^ 

k(^.2) ^ 0{e~^^'' ), which can be neglected in comparison 
with . 

B. Consider now the case where the dimension- 
less intergrain tunneling conductance g is large, g = 
7r\t\'^i/{eF)'^ ■ One can show 24] that the diamagnetic 
contribution to the electrical conductivity is cr^''-' = 

e^a^-'^g (^1 - ^ In . Comparing with Eq.||ni), it is 
evident that for the diamagnetic components, the kinetic 
part of the thermal conductivity is simply related to the 

jS", where 



electrical conductivity by k'^'^'-^^ /(Ta^'^^) 
we have restored the Boltzmann constant ks- So at this 
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level the WF law is obeyed. Next consider the kinetic 
part of the paramagnetic contribution to the thermal 
conductivity, k(P'-^)(T). Ref.^ shows that cr(P) < 
adopting a similar argument, one can show that k^P'^'> 
is smaller than n^'^'^') by a factor g~'^, and furthermore 
it does not have a logarithmic singularity. We have not 
yet considered the extra contributions from the potential 
part, K^'^'^^ and k'^P'^\ Evaluating k^'''^' for example 
leads to corrections similar to those obtained in Ref.[l3|- 
However since the potential part k^^'^^ of the paramag- 
netic term is of the order of k'-'^'^-', a proper treatment 
should consider both contributions. This could, in prin- 
ciple, lead to a different conclusion from Ref.[l^. So at 
this stage we are only able to make a weaker statement 
that for the strong tunneling case at not too low tem- 
peratures {T > gE^e-^az^, ^ = ^ + 0{g-^). Thus 
the Lorenz number L is larger than Sommerfeld's free 
electron value by 0{g~^) corrections. 

We conclude with a discussion of our results and com- 
ments on existing and future experiment. We have 
shown that thermal transport in granular metals at 
low temperatures is dominated by cotunneling of low- 
energy electron-hole pairs. These neutral excitations do 
not transport charge. As a result, while the electri- 
cal conductivity for weakly coupled grains is exponen- 
tially small in temperature, cr(T) ~ e^ga^~'^e~^''^^ , the 
thermal conductivity is only algebraically small, k{T) ~ 
a'^~'^k^g^T^ / E^. The particle-hole cotunneling process 
is physically equivalent to a particle cotunneling loop. 
An electron executing a cotunneling loop brings back 
its charge to the starting grain and hence there is little 
change in the electrical conductivity. There is no require- 
ment, however, that the returning electron has exactly 
the same energy. This conservation of grain charge but 
not grain energy in a cotunneling loop is at the heart of 
the difference between heat and charge transport. The 
energy gained by an electron due to nearest neighbor tun- 
nchng is of the order of (t'^/Ec), and the tunnehng prob- 
ability is, roughly, [t/EcY. The number of particle-hole 
excitations in each grain is of the order of i^{ep)T ^ 1, 
and the energy of a particle-hole is of the order of T. This 
gives us a rough estimate k(T) oc g'^T^/E'^, in agreement 
with our detailed calculation. The near-continuum of 
states in the metallic grains means the electron-hole ex- 
citations can have very low energy; in this sense cotun- 
neling in regular granular metals differs fundamentally 
from exciton transport in the Hubbard model even when 
intergrain tunneling is small. However there seems to be 
some resonance with excess heat transport by delocalized 
tightly-bound excitons in dirty semiconductors 0, . 
Experimentally, Al — Ge granular systems in the g S> 1 
regime show a linear— T thermal conductivity with 
the value of the Lorenz number somewhat larger than Lq, 
L « 1.81/0 at the lowest temperatures measured. This is 
consistent with our assessment. 



Our predictions can be experimentally verified: even 
though the g^T^/E'^ law we obtain for k(T) is remi- 
niscent of the phonon contribution, it can be dis- 
tinguished experimentally through its dependence on g 
and Ec- Besides, the electronic contribution may be of 
the order of the phonon contribution as can be seen in 
the following rough estimate. Let a denote the distance 
between grain centers, I the phonon mean free path, c 
the speed of sound, n the atomic number density, and 
Qd the Debye temperature of the material. The ratio 
of the electronic and phononic thermal conductivities is 
then of the order of r = {gQD/Ec)'^{kBQD)/{fi(incl). In 
a typical three-dimensional granular metal, the grain di- 
ameter is of the order of 100 A and the insulating space 
between the grains is about lOA. With Qd 300K, 
Ec ~ lOOK, 5 < 1, a - lOOA, I - lOAjl^, c ~ lO^m/s, 
and n ~ lO^^m^'^, we estimate r k g'^ ^ 0(1). The 
electronic contribution can be further enhanced by in- 
creasing the size of the grains and choosing an insulator 
with a higher dielectric constant; each of these decreases 
Ec- Physically, the surprisingly large electronic thermal 
conductivity can be attributed to the large electron speed 
compared to the speed of sound. 
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